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Abstract. - We study disorder-induced spectral correlations and their effect on the magnetic 
susceptibility of mesoscopic quantum systems in the non-diffusive regime. By combining a dia- 
grammatic perturbative approach with semiclassical techniques we perform impurity averaging 
for non-translational invariant systems. This allows us to study the crossover from clean to 
diffusive systems. As an application we consider the susceptibility of non-interacting electrons 
in a ballistic microstructure in the presence of weak disorder. We present numerical results for 
a square billiard and approximate analytic results for generic chaotic geometries. We show that 
for the elastic mean free path I larger than the system size, there are two distinct regimes of 
behaviour depending on the relative magnitudes of I and an inelastic scattering length. 



Phase-coherent, disordered conductors where the electron motion equals a random walk 
between impurities have traditionally been of interest in mesoscopic physics jl]] . Random walks 
occur in the diffusive regime where the elastic mean free path £ is much smaller than the system 
size L. On the other hand, the development of high-mobility semiconductor heterostructures, 
combined with advanced lithographic techniques, have allowed the confinement of electrons 
to two-dimensional microstructures of controllable, non-random geometry. They have been 
coined "ballistic" since £ > L. Nevertheless, residual impurity scattering is nearly unavoidable 
even in these systems and it became clear that disorder can be strong enough to mix energy 
levels and effect the two-level correlation function, K(ei,£2), even in the "ballistic" regime 

p. There it is necessary to consider both disorder averaging, (. . .) d , and size (or energy) 
averaging, {. . .) L . After such averaging the two-level correlation function may be divided 
into two separate terms [§, ( K d (e 1 , e 2 ) ) L = {{v{ei)v{e2))d) l - {{v{ei))d{v{e 2 ))d) l and 
K l (e\,E2) — {{v{E\))d{v{£2))d)L — v 1 , where v denotes the single particle density of states 
and v = ((i/(e)}d} l its mean part. K d is a measure of disorder-induced correlations of v, while 
K L is given by size-induced correlations. 

The orbital magnetism of isolated mesoscopic systems has been the subject of much theo- 
retical interest, in particular because it is sensitive to spectral correlations: For a system with 
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a fixed number of particles it is necessary to consider averaging under canonical conditions 
@> H> D resulting in a large contribution to the average magnetism. The corresponding sus- 
ceptibility is given by |§] (x(H)) = - (A/2) d 2 /dH 2 (SN 2 (fi; H)) . Here H is the magnetic 
field, A is the mean level spacing and (5 N 2 (/x; H)\ is the variance in the number of energy 
levels within an energy interval of width equal to the chemical potential, /i. This variance is 
related to K(si,E2',H) by integration of the level energies s\, £2 over the energy interval. In 
the following we will label the contributions to the susceptibility, corresponding to ^ K d V and 
K L , as ( X d (H)) and ( X L (H)), respectively, so that ( X (H)) = ( X d (H)) + ( X L (H)) ■ 

For real systems, besides £, there are additional relevant lengthscales at which inelastic 
scattering (L^) or temperature smearing {Lt) produce a damping of propagation. For clarity 
we refer to such a lengthscale as L^, although we assume that similar general arguments will 
hold for finite Lt- For ballistic motion, L$ is related to the level broadening 7 by 



L<j, kpL A 



(1) 



L 2?r 7 

with Icf being the Fermi momentum. It divides the "ballistic" regime into two sub-regimes. 
In the first, L, L^ < £, the particle motion is nearly ballistic since damping due to inelastic 
scattering typically occurs before impurity scattering; for the remainder of this paper we refer 
to this regime as inelastic. In the second, L < £ < L^, a particle may scatter many times off 
impurities before scattering inclastically and we refer to this regime as elastic. 

In this paper we employ a systematic semiclassical approach to calculate the contribution of 
orbits of all lengths in the elastic and inelastic regimes to weak field (x d {H)) for microstruc- 
tures with white noise disorder 0. Energy and impurity averaging are performed within 
a diagrammatic perturbation approach applicable to the "ballistic" regime. In contrast to 
many techniques valid for bulk disordered systems, £ < L, we do not assume translational 
invariance. Instead, we use an approach related to the "method of trajectories" devised for 
thin superconducting films H and we write Green functions semiclassically in terms of classical 
paths which include the effect of boundary scattering. Note that similar |l, and alternative 
[ [To| , methods have been used to consider weak localization in thin films in a parallel magnetic 
field. Our approach allows us to study the complete crossover from diffusive to clean systems 
for arbitrary values of (smaller than Vptu, where tn is the Heisenberg time). We present 
results for ballistic systems with both chaotic and integrable dynamics in the clean limit. As an 
example of an integrable geometry, we treat the case of the square billiard. Experiments |TT| 
on the orbital magnetism of ensembles of squares were performed in the "ballistic" (inelastic) 
regime, motivating theoretical studies of the susceptibility for L < £ jl2|, |l3| . 

Semiclassical diagrammatic approach. - We begin by presenting some more details concern- 
ing our semiclassical evaluation of the disorder correlation function, K d , and corresponding sus- 
ceptibility ( x x d {H)'). We consider non- interacting electrons in a weak, perpendicular magnetic 
field. In terms of retarded and advanced single particle Green functions, C? + ( - ' (ri, r 2 ; e; H), 
K d may be written as K d ( £l , e 2 \ H) « (A 2 /2tt 2 ) K ((tr H)tr G~{e 2 ; H))) d , where ((. . .)) d 

implies the inclusion of connected diagrams only. Using a diagrammatic approach introduced 
by Altland and Gefen ||] the field sensitive part of K d can be expressed as a sum including 
Cooperon type diagrams defined by 

r ~\rb I 

SP(lo;H) = Tr [c (c >] = / Y[d d r, \{ ^ C \r m ,r m+1 -cj;H) ; r n+1 = n . (2) 

J j=l m=l 

Here C (C) (ri, r 2 ; u; H) = (1/2ttPt) G+(n, r 2 ; e i; fl)G-(n, r 2 ; e 2 ; H), = (G +( -- ) ) d is the 

disorder averaged single particle Green function, u> — e\ — £2, and r = £jvp. 
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Semiclassically, G + (ri,r 2 ) can be expressed as a sum over classical trajectories t between 
ri and r 2 O, 



G + (r!,r 2 ) ~ ^D t (ri,r 2 )exp 



1 q< r r ) £t(ri ' r2) 



(3) 



The prefactor D t includes the classical phase space density, St stands for the classical action 
along an orbit t (in the absence of disorder) including the Maslov index, and L t is the orbit 
length. C r 2! ^\ H) is then given in terms of pairs of classical paths which explicitly 

include the effect of boundary scattering. However most pairs (of different paths) produce 
oscillating contributions which we assume to vanish after energy or size averaging ( 1 ). The 



main contribution to the held sensitive part of ( K d J L arises from diagonal terms (otherwise 
known as the Cooperon channel) obtained by pairing paths with their time reverse. Assuming 
that the magnetic held affects the phase of the particles but not their trajectories we can write 
C^(ri,r 2 ;o;;fO = Eta™ C,¥\r u r 2] u>;H) where 



v F \Dt\ 

2-KVi 



■ exp 



La 



Lt 

I 



+ iuT* + i- 



.47T 

Vo 



A.dr 



(4) 



Here, T t is the period of the trajectory, A is the vector potential, tpo — hc/e, and the level 
broadening was introduced via to — > u> + 17. Eq. (Q) depends, besides t, only on the system 
without disorder and holds for both integrable and chaotic geometries. 

The disorder induced contribution to the average susceptibility is given by ((p = HL 2 /tpo): 



XL 



6 d 2 

TT 2 dip 2 



00 -. 



(5) 



71=1 



where the bulk Landau susceptibility is \l — — e 2 /247rmc 2 for spinless electrons. In Eq. (|^) 
the S { n ] are now assumed to contain diagonal terms only and their contributions (Eq. (||)) can 
be calculated by diagonalising C,^ which in general cannot be done analytically. However we 
can use the fact that all the variations of C}- c ^ occur on classical lcngthscales; rapid oscillations 
on the scale of A_f cancel out. It is thus possible to discretise the "classical" operator on a 
lattice in space with grid size greater than Af- By summing over trajectories betwe en l attice 
cells one can compute C^ C \ an d thereby (x d (-ff)), efficiently by numerical means [TilP^I . 

The propagator (above Eq. (|J)) is made up of a summation over all diagonal pairs 
of paths (including boundary scattering) between any two given impurities situated at ri 
and r 2 . On taking the trace over n propagators C , one sees that the field sensitive part 
of S n , Eq. (|^), consists of a summation over flux-enclosing closed pairs of paths (in position 
space) involving n impurities and an arbitrary number of boundary scattering events. However 
this summation does not include closed pairs of paths which follow periodic orbits of the 
corresponding clean system. Such paths involve zero momentum transfer between the Green 
functions at the impurity positions; they actually represent disconnected diagrams which are 
included in {y^iH)) and must not be counted again in the determination of (x d {H))- It is 
the presence of these periodic orbits in the determination of (x L {H)) (_] that leads to strong 
sensitivity with respect to the system geometry [|l(| (see below). 

( 1 ) In a ballistic system additional oscillatory terms will however remain upon pure disorder average 
for fixed size HQ. 

( 2 ) A similar calculation has been applied to a different "classical" operator, describing interaction 
effects in ballistic quantum dots, in Ref. ]l5| . 

( 3 ) A similar sensitivity with respect to the system geometr y o ccurs in interacting ballistic quantum 
dots due to the presence of off-diagonal periodic orbit terms 
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Fig. 1. - Disorder-induced average susceptibility (x d (0)) for a square geometry in the elastic regime 
(L < £ < L$) as a function of the elastic mean free path I for UfL — 60 and two strengths of inelastic 
scattering, 7/A = 1 (lower), which corresponds to L^/L » 9.5, and 7/A = 0.392 (upper). The dashed 
horizontal line indicates the result by GBM for 7/ A = 1. The inset shows (x d (0)) as a function of 
UfL for 7/A = 1. From the top, the 5 curves are for values of i/L = 2, 4, 5, 1 and 10. 

Fig. 2. - (x d (0)) in the inelastic regime {L,L^ < i) as a function of I for kpL — 60 and 7/A = 10 
which corresponds to L^/L ~ 0.95. Circles are our numerical results and the triangles are numerical 
results including only the contribution of Si. The solid and dashed lines are both fits (see main text). 
The inset shows (x d (0)) as a function of UfL for 7/A = 10. From the top, the curves are for values 
of i/L = 2, 4 and 8. The symbols correspond to our numerical results and the solid lines to a fit. 



In the following we apply the above formalism to the case of an ensemble of disordered 
square billiards. Gefen, Braun, and Montambaux (GBM) considered the contribution 
of trajectories longer than I to (x d (-ff)) in an approximate way, while Richter, Ullmo, and 
Jalabert (RUJ) jl3| calculated (x L (H)) for a square by assuming that the disorder perturbs 
the phase, but not the trajectory, of semiclassical paths of the corresponding clean geometry. 
We first compute the (Eq. (||)) for the square geometry by employing the extended zone 
scheme fu| to write £^ C K r ij r 2 ; u>; H) as a sum of propagators along straight line paths. We 
then perform a complete calculation of (x d (H)), Eq. (||), in the elastic and inelastic regimes 
and compare the results with those by GBM and RUJ. 

Elastic regime: L < £ < V - Fig. 1 shows (x d (0)) as a function of £ for a typical 
experimental value of fc^ L = 60. The lower curve is for 7/A = 1 (i.e. L^/L w 9.5 at UfL = 60) 
and the upper is for 7/A = 0.392 p8| . Note that in the diffusive regime, £ < L, there is linear 
increase with I in agreement with Ref. ]l9[ ]. For L < £ < L^, we find a weak dependence 
of (x d (0)) on i. Our result is on the whole close to the prediction by GBM [|l2| who found 
a paramagnetic ^-independent contribution, (x d (0)) / \xl\ ~ 0.23 kpL(A/j), shown as the 
dashed horizontal line in Fig. 1 for 7/A = 1. As I increases further, (x d (0)) decreases; we 
discuss this behaviour in more detail later when considering the inelastic regime. 

The inset of Fig. 1 shows (x d (0)) as a function of UfL for 7/A = 1. From the top, the 
five curves are for values of i/L = 2, 4, 5, 1 and 10. For all disorder strengths, (x d (0)) is 
paramagnetic and it increases linearly with UfL. For L < £ < L^ the gradient of the curves 
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Fig. 3. - Comparison of (x d (0)) and (x L (0)) for the square. The solid line shows our numerical 
results for (x d (0)) and the dashed line is an analytical expression for (x L (0)) (see main text) as a 
function of L^/L for UfL = 60 and £/L = 2. Inset: value of L^/L at which the two contributions 
are equal as a function of £/L. The solid line is the analytical estimate Eq. (^) and the circles are 
obtained by comparing our numerical results for (x d (0)) with the analytic expression for \X L (0))- 



Fig. 4. - Comparison of the (normalized) semiclassical estimates for (x d (0)) (solid line) and (x L (0)) 
(dashed) for a generic chaotic geometry for UfL = 60 (see main text). Inset: the (straight) line shows 
the values in the (£, L^)-plane where both contributions are equal. 



is approximately independent of I and we find it to be ~ O.I8A/7. However there is a fcpL 
independent offset to the curves which is £ dependent and not described by GBM. 

Inelastic regime: L,L4, < £. - Fig. 2 shows (x d (0)) as a function of £ for fcpL = 60 
and 7/A = 10 (i.e. L^/L sa 0.95). Circular points correspond to our full numerical results, 
while triangular points represent the contribution of S\ only (minus the disconnected part) 
in Eq. (||). The solid line is the equation (x d (0)) = (a>oL/£) exp (— aiL/£) where ciq and 
ai are fitting parameters. The dashed line displays an equation of the same type but with 
a\ = 2\/2 and ao as the only fitting parameter. For £ 3> L the susceptibility is dominated 
by the contribution with the lowest number of impurity scatterings S±. As £/L is reduced, 
progressively more terms in the summation of Eq. (||) become relevant, and there is good 
agreement with the solid line fit for £ > L. 

The inset of Fig. 2 shows (x d (0)) as a function of UfL for 7/A = 10. From the top, the 
three curves are for values of l/L = 2, 4 and 8. The symbols correspond to our results and 
the solid lines to (x d (0)) = 60 exp [—bi(j/A)/kpL] where 60 and b\ are fitting parameters. 

Comparison with the contribution of clean correlations. We compare the magnitude 
of <x d (0)) with that of <x L (0)) for squares. It has been shown |l^, |l^] that the low field 
susceptibility of an ensemble of clean squares is dominated by the shortest flux enclosing 
periodic orbits of length L t = 2^/2L and their repetitions over a broad range of temperature 
(and thus inelastic scattering strengths). For the ballistic white noise case considered here, 
the effect of disorder averaging on the susceptibility was described by an additional damping 
exp(— L t /£) of the response of the clean system This result corresponds to (x L (0)) 

including the disorder damping exp(— L t /2£) of the single particle Green functions. 
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We use the results of RUJ |13|, at zero temperature and introduce the level broadening 
7 in the same way as for the disorder correlations above |^0[ . It is then possible to sum the 
contribution of all repetitions of the fundamental orbit explicitly which gives (x L (0)) / \xl\ — 
(V2/bTr)k F L/ sinh 2 [V2(L/e + L/L^)}. Fig. 3 shows this expression for (x L (0)) (dashed line) 
and our numerical results for (x d (0)) (solid line) as a function of L^/L for kpL — 60 and 
i/L = 2. Although (x d {0)) < (x L (0)) I0r small L^/L and vice versa for large L^/L, it is clear 
that both contributions are relevant over a broad range of L$/L. We make an estimate for 
the value of L^/L at which the contributions of (x L (0)) and (x d (0)) are equal by comparing 
the above analytic approximation with that given by GBM. We find for I > L that (x d {H)) 
is larger than (x L (-ff)) f° r ^<f> greater than a crossover value, IA, given by 



L% {k F Lf 



L 2tt 2 



v2,v/ £ ( ~) +8ir 2 (~ 



(6) 



The inset of Fig. 3 shows this estimate (solid line) compared to points (circles) obtained by 
comparing our numerical results for (x d (0)) with the analytic expression for (x L (0))- 

The experiment on the orbital magnetism of ensembles of squares had estimated values 
for the elastic mean free path of i/L ~ 1 — 2, for the phase-coherence length of ~ (3 — 10)L and 
for the thermal cutoff length of Lt/L ~ 2. Hence, the lengthscale L^ (Eq. ([|)) is determined 
by the shorter length Lt- Fig. 3 shows that for these experimental parameters (and for white 
noise disorder) both the disorder and size- induced correlations are relevant, however the latter 
contribution is dominant. The measured value of the susceptibility at low temperature was 
x(0) ~ 100|xl|, with an uncertainty of about a factor of four. After including a spin factor 
of 2, the combined contributions (x d ) and (x L ) calculated above, together with an interaction 
contribution of the same order [jlB), are in broad agreement with the experimental result. We 
note, however, that a theoretical explanation of the temperature dependence of the measured 
susceptibility is still lacking. 



Experimental ballistic structures as those in Ref. 1 1 are usually characterised by smooth 
disorder potentials. The effect of smooth disorder on (\ L ) has been analysed in Ref. 1 1 3 
showing that the reduction of the clean contribution is less strong than for white noise disorder 
and no longer exponential. Smooth disorder effects can be incorporated into the present 
calculation by introducing an angle-dependent cross section for the impurity scattering between 
two trajectory segments. 

Chaotic geometries. - For systems with a generic chaotic, clean counterpart we obtain an 
analytical estimate for (x d (0)) in the elastic regime after transforming the sum over densities 
\D t \ 2 in Eq. ([|) into probabilities P(r, r';i|A) to propagate classically from r to r' at time t 
accumulating an "area" A [Q. Assuming a Gaussian "area" distribution with a variance a, 
which is taken to be t independent, we find 

5r)( W ;F)«{l + ^^ + ( 7 -f w) r}" n . (7) 



Summation of the <S„ , Eq. (||), leads to 



\ X l\ L\L^ + i) 



(8) 



(x d (0)) is shown as the full line in Fig. 4 which is rather similar to the corresponding curve 
in Fig. 3 for the square geometry. A corresponding approximation for (x L (0)) |17| is shown 
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as the dashed line in Fig. 4. Both can be shown to add up to an ^-independent response 
(x(0)) /\xl\ — 96(rL l p/L i for generic chaotic geometries. 

Conclusion. Disorder-induced spectral correlations and their effect on the magnetic 
susceptibility in the non-diffusive regime £ > L were considered. We focused on the square 
billiard, whose corresponding clean geometry is integrable, and showed that there are two 
distinct regimes of behaviour depending on the relative magnitudes of I and L^. This approach 
enabled us to study the complete crossover from diffusive to clean systems for arbitrary values 
of i and (smaller than Vptjf). Note that it may be possible to calculate the susceptibility for 
values of La and I greater than vptfj using a non-perturbative approach such as the ballistic 
a model 
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